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An open quantum system interacting with its environment can be modeled under suitable as- 
sumptions as a Markov process, described by a Lindblad master equation. In this work, we derive 
a general set of fluctuation relations for systems governed by a Lindblad equation. These identities 
provide quantum versions of Jarzynski-Hatano-Sasa and Crooks relations. In the linear response 
regime, these fluctuation relations yield a fluctuation-dissipation theorem (FDT) valid for a station- 
ary state arbitrarily far from equilibrium. For a closed system, this FDT reduces to the celebrated 
Callen-Welton-Kubo formula. 

PACS numbers: 



I. INTRODUCTION 

Fluctuations in non-equilibrium systems have been shown to satisfy various remarkable relations [2(3, [2l|, [3(| [4l|, |47], 
EH H3] discovered during the last twenty years. These results have lead to fierce discussions concerning the nature 
of heat, work and entropy, raising the fundamental issue of understanding the interactions between a given system 
and its environment (e.g., a thermal bath). In the classical realm, these problems have been progressively clarified 
whereas they arc still under investigation in the quantum world. 

Historically, quantum fluctuation relations were first studied by Callen and Welton in 1950 [l2j|. These authors 
derived a Fluctuation-Dissipation Theorem for a closed quantum system isolated from its environment, initially 
in thermal equilibrium and suddenly perturbed by a small time-dependent term added to the time-independent 
Hamiltonian H. This approach was further systematized by Kubo [5l| in the linear response theory. 

Since 2000, three main directions have emerged for constructing quantum fluctuation relations which extend the 
linear response theory of fl2l . [5l| and which are quantum analogues to fluctuation relations for classical systems. In 
the following, we rapidly review these three different routes by emphasizing their goals and differences. We also cite 
useful articles. 

1. In the first approach, initiated in 2000 by Yukawa [76| and continued by Allahverdyan and Nieuwenhuizen Q, 
a definition for the quantum work operator is introduced but this operator does not obey any fluctuation relation. 
More precisely, for closed but non-autonomous systems (with time-dependent Hamiltonian H t ), it is proved in 0] 
that the naive (inclusive) work operator Hf 1 (t) — (0) (where the exponent H stands for Heisenberg picture) does 
not satisfy any fluctuation relation, unless [Ho,H t ] = 0. This implies the existence of quantum corrections to the 
Jarzynski relation fvj\ . The works of Bochkov and Kuzovlev in 1977 [l(|, and of Stratonovich in 1994 (67} can be 
viewed as earlier attempts in this direction. 

2. The second line of attack was opened by Kurchan and Tasaki in 2000 j53l[73j. and was continued by many groups 
0, EH HH IH, l70l - [72T ] . This approach is based on a different definition for the work. The idea is to introduce initial 
and final measurements of the system's energy, according to the quantum-mechanical measurement postulates. The 
work (which is viewed as an energy difference) is then a two-point functional obtained by calculating the difference 
between the results of the two measurements. This definition of the work differs fundamentally from the previous 
one (see [26l [33j for a comparison) and does satisfy quantum fluctuation relations. The results of this approach are 
thoroughly reviewed in [l6|, HH . In summary, we can say that this definition of the work as a two-point measurement 
has been applied in the following different contexts: 

(i) Closed but non-autonomous systems prepared in a Gibbs state and isolated from their environment during their 
evolution (which is thus unitary). Kurchan [53[ studied the time cyclic case with H t = Hx-t (T being the period) 
and proved Jarz yns ki and Crooks relations in this set-up. Tasaki [73j generalized this result to the non-cyclic case. 
Many groups fj[, l6ll. |62|. |7(I l7l| simplified the theory. In particular, in 2007, Talkner et al. [7(3] clarified the fact that 
work, being characterized through a process with two measurements, 'is not an observable' and cannot be represented 
by any hermitian operator. The link and the difference with the older Bochkov-Kuzovlev approach [Tol | (inclusive 
versus exclusive work) was explained by Campisi et al. in I15| and the possibility to perturb the unitary evolution 
by N-points measurements was studied by Callens et al. in [13j . Finally, a more general relation for a 'quantum 
generating functional' was derived by Andrieux et al. (see equation 12 of Q)- 
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(ii) The second set-up corresponds to the general case of an open system continuously interacting with its surround- 
ings. This case can be formally reduced to the previous situation considering the system together with its environment 
to be a global, closed, system [l4|, [H, [H, [H, [30, [H, [6l|, Iz3 • The main physical advantage of this approach is that one 
can use the previous expression of work as a 2-points measurements. But, such a 'holistic' approach, involving both 
the system and its surroundings, leads to fluctuation relations that are difficult to assess experimentally. Indeed, the 
environment is usually large and hardly controllable and only the degrees of freedom are experimentally accessible. 

(iii) Finally, the study of heat-matter exchange for two systems in contact Hfl Hll H^, [65], [66| can be viewed as a 
special instance of the previous set-up. Contact between reservoirs at different temperatures and chemical potentials 
lead to transport of energy and matter. A famous example is electron counting statistics [50j in which small nanoscale 
electronic devices exchange electrons. 

3. In a third category of works, the system is modeled by an effective master equation at the mesoscopic scale. 
This approach was pursued in relatively few articles [22|, [28|, HH [76| . The idea is to consider an open system 
continuously interacting with its environment and to project out the degrees of freedom of the bath to derive an 
effective dynamics for the system. Then, it is assumed that the reduced dynamics is described by a closed evolution 
equation for the density matrix of the system. Under some further assumption s l42l. l56l , this master equation can be 
brought into a Markovian form known as the Lindblad equation. In the papers |30l.l34|. the quantum master equation 
is treated as an effective classical master equation associated to a pure-jump process; this allows the authors to use 
the concept of pathwise trajectory and the trajectorial definition of entropy production. The classical fluctuation 
relations can then be applied. This approach [30, [H|, albeit very powerful because it rests on the highly developed 
field of classical fluctuation relations, doesn't provide us with any explicit relation involving quantum observables. 
Quantities such as entropy or work are defined along effective classical trajectories and their transposition for the 
original quantum system are not at all obvious (and in fact, the effective classical process is not uniquely defined). 
A different philosophy, adopted in [22T [28l |76| . is to work directly with th e q uantum master equation. However, in 
[22T |76| the time evolution was discretized in an ad-hoc manner, and in [28| the transition rates were given in an 
arbitrary way. At a conceptual level, one could object the relevance of the approach 3, valid only at a mesoscopic 
scale. Since fluctuation relations can be established at the microscopic scale (by the approach 2), the theory at 
the mesoscopic scale should simply result from the proved microscopic relations. This objection is not valid for two 
reasons. First, one should recall that even in the classical case, the fluctuation relations were experimentally tested 
for effective stochastic models, valid only at a mesoscopic scale A second argument, given by De Roeck in 

[lo], emphasizes the fact that for a mesoscopic quantum system entropy production is not well defined (in contrast 
to the case of a classical system) mesoscopic fluctuation relation can not be obtained by a coarse-graining procedure. 



The present work follows the third approach to quantum fluctuation relations. We study the non-equilibrium 
fluctuations of an effective open quantum system modeled by a Lindblad master equation. The Lindbladian evolution 
is a non-unitary dynamics for the density matrix p t of the system, described by a differential equation with a generator 
L t (semi-group property). This effective Markovian description is widely used in Quantum Optics [42jj. But unlike 
[22I l28l . [30, [3J, [zol, our goal, here, is to work directly with the continuous time Lindblad equation, to define an 
associated time-reversed dynamics and to derive fluctuation relations with quantum observables. Therefore the 
fluctuations relations we obtain stem from structural and symmetry properties of the Lindblad master equation. 
The key results of the present work, given in Eqs. (|50I5II22I39|) . represent an original contribution to quantum non- 
equilibrium statistical mechanics. 

Our strategy will be to use a suitable deformation of the master equation, which will allow us to prove a generic 
relation amongst correlation functions, a kind of book-keeping formula which is a quantum analog of Jarzynski-Hatano- 
Sasa relation and Crooks relation. Furthermore, by a lowest order expansion, we derive a generalized fluctuation- 
dissipation theorem valid in the vicinity of a quantum non-equilibrium steady state. For the special case of a closed 
system, our approach retrieves previously known work identities @, HH, [5^, [70, [zl, [Zll as well as the quantum equi- 
librium fluctuation-dissipation theorem [l2|, [HI ■ 

The outline of this work is as follows. In section [ill we recall some basic properties of the Lindblad equation 
which represents an effective Markov evolution for a system in contact with an environment. We write a formal 
solution for the Lindblad equation and give the expression of multi-time correlation functions. In section Hill we prove 
the Quantum Jarzynski-Hatano-Sasa relation associated with the Lindblad equation. We use this relation to derive 
a steady-state quantum Fluctuation-Dissipation theorem that generalizes the Kubo-Callen-Welton formula (which 
applies only to a closed system in the vicinity of thermal equilibrium). In section ITVl we investigate the properties 
of a given Lindblad evolution under time-reversal. This allows us to prove for open quantum systems a general 
version of the Tasaki-Crooks relation. Concluding remarks are given in section [Vl Technical details are deferred to 
the appendices. 
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II. MASTER EQUATION FOR QUANTUM MARKOV DYNAMICS 

Consider a quantum system S in contact with a thermal reservoir (or environment) R. The total system S + R 
is closed. Its evolution is unitary and the total density matrix evolves according to the quantum Liouville equation 
governed by the total Hamiltonian -fftotai which can be broken into three pieces: the Hamiltonian of the (small) system 
Hs, the Hamiltonian of the reservoir Hr and the interaction Hamiltonian between the system and its environment 
Hi. We are interested in the degrees of freedom of S and therefore we would like to have at our disposal an evolution 
equation for the density matrix p t of the quantum system S alone, the degrees of freedom of the reservoir being traced 
out. Generically, such an equation is non-local in time: the coupling Hi induces memory effects. However, under 
some specific assumptions, a differential equation of first order with respect to time can be derived for p t : one must 
assume that (i) the full system S + R is prepared in a correlation free state; (ii) the reservoir R is large enough so 
that it has a very short memory time r c (more precisely one must have r c -C h/\Hi\ where \Hi\ is the typical order 



of magnitude of the interaction Hi matrix elements 4J|). This second assumption is the crucial Markov hypothesis: 
when it is fulfilled a memoryless and coarse-grained description of the system S becomes possible. This condition is 
generally satisfied in Quantum Optics [U, |44| . 

There are two main methods to derive the master equation for quantum Markovian dynamics. One way is to 
make a precise model for the reservoir (typically an infinite set of quantum oscillators) and to eliminate explicitly the 
environmental degrees of freedom. The Markov approximation can be analyzed and justified precisely [2J]. Another 
possibility is to study the structural properties of the 'quantum map' that carries the density matrix p t at time t 
into the density matrix pt+dt at time t + dt. Such a map must be linear, hermiticity preserving, trace conserving and 
positive. In fact, this map lifted to an operator on the total system S + R (for any given environment) must remain 
positive (this stronger requirement is called complete positivity) p2l . l44j . These physically reasonable conditions arc 
stringent enough to mathematically constrain the possible forms of the evolution equation of a quantum Markov 
dynamical system (5(| . The resulting equation is called a Lindblad equation. Its generic form, equation ([1} , and some 
of its basic properties are discussed below. 

A. Some properties of the Lindblad equation 

We consider a quantum system prepared initially with a density matrix ttq. Because of its interactions with its 
environment, the density matrix of the system, becomes a function of time, and will be denoted by p t . In the present 
work, we adopt the framework of Quantum Markovian Dynamics. The evolution of pt is thus given by a Lindblad 
master equation, which can be written in the generic form [42l |44| : 



d tPt = -i[H u Pt] + J2 ( V ^ V i - lVi V iPt ~ \ptV}V?) . (1) 



On the right-hand side of this equation, the first term —i[H t , X] is the conservative part where H t is the Hamiltonian of 
the system that may depend on time. The other terms represent the interactions of the system with its environment 
(also called the 'bath') and also represent the effect of measurements (i.e. dissipation and (de)coherence effects). 
The operators Vi are called the Kraus operators, they are not necessarily hermitians and may depend explicitly on 
time. The Kraus number I depends on the system considered. If the system under consideration is closed, all Kraus 
operators vanish identically, Vi = 0, and the Lindblad master equation reduces to the quantum version of the Liouville 
equation. Equation ([T]) can be written symbolically as 

dtPt = L\(H , (2) 

where we have introduced the Lindbladian superoperator L\ which acts on the density matrix p t and generates its 
time-dynamics. We emphasize that L\ is a superoperator because it is a linear map in the space of operators. The 
fact that we have used in equation @ the symbol L t for the Lindbladian rather than the more usual notation 
L t is purely a matter of convention: this will allow us to write some expressions of time-ordered correlations of 
observables in a simpler manner (because L\ acts on the density matrix p t and its conjugate superoperator L t acts 
on observables) . More precisely, the space of operators is endowed with the following Hilbert-Schmidt scalar product 
(Y,X) = Tr(Y'X), where X and Y are arbitrary operators and Y* is the hermitian conjugate of Y. This allows us 
to define a pair of adjoint superoperators L t and L t as follows 

(Y,L t X) = Tr(Y\L t X)) = (l\y,X) = Tr{{L\Y)^X) . (3) 
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A simple calculation allows us to write 

LtX - i[H t ,X] + (vJXVi - \v}ViX - \XV?V?) . (4) 
i=l ^ ' 

Two important properties of the Lindbladian L t are L t l = (Trace conservation) and L t (X^) = (L t Xy. 

The Lindblad equation is extensively used in Quantum Optics. A simple example is a two-level atom emitting a 
photon in free space. The density matrix p t is a 2 by 2 matrix and the Kraus operators reduce to Pauli lowering and 
rising operators. The Lindblad equation is then simply a set of four coupled first order differential equations |44l |. 



B. A formal solution of the Lindblad equation 

The quantum Master equation can be solved formally by introducing the evolution superoperator Pq: 

Pt={P&n , (5) 

where ttq represents the density-matrix at the initial time. The evolution superoperator between the two times 
s < t is defined by 

P* = / duL u =1 + V / JdU L tl L t2 ...L tn . (6) 

V-^s / „ =1 •/s<ti<t 2 <...<t„<t i=1 

In this time-ordered exponential, time is increasing from left to right. This symbolic writing will be very useful to 
perform formal calculations and to write perturbative expansions. Let us recall how equation ([3]) is proved. First, 
we observe that this equation is true at t = because P® is the identity operator. Then, from the time-ordered 
exponential ([5]), we find f t Pl = P l s L t . This leads us finally to 

j t p t =(L\{P^)^=L\ Pt . (7) 

Thus, pt satisfies the Lindblad equation @ with initial condition tto. We note that this technique of proving an 
identity between operators, such as equation ([5]), by showing that both operators are solutions of the same (first 
order) differential equation with the same initial condition, will be used repeatedly in this work. 



C. Expression for multi-time correlations 



Using the evolution operator, we can write a general expression for multi-time correlations of different observables. 
For < ti < t 2 < ■ ■ ■ < ijv < t, the time-ordered correlation of observables Oq,Oi,C>2---On is given by 

(Oi(ii)0 2 (i 2 ) ■ • • N (t N )) no = Tr (noP^O^Oi . . . P^On) ■ (8) 

A justification of this expression can be found in [42| or in [IT| : One starts by the correlation of two operators at two 
different times. This correlation can be evaluated in the Heisenberg representation by using the full Hamiltonian of 
the system plus its environment. In order to obtain an expression that refers to the system alone, the partial trace 
over the environment has to be performed. Using the same assumptions (factorization of the initial condition and 
weak coupling) that lead to the Lindblad Master equation ©, one shows that the time-ordered two-time correlation 
function satisfies an evolution equation which is the dual of the equation ([2]), i.e. it is governed by the Lindblad 
operator L t . This proves the formula for N = 2. The general case is then obtained by induction. Besides, a 
rigorous proof of the formula (O is given in for a toy model of one harmonic oscillator in a bath of independent 
harmonic oscillators. 

We recall that, in the formula the operator ttq represents the initial density matrix of the system and emphasize 
that the superoperator p** +1 operates on all the terms to its right. For example, for N = 3, the explicit expression, 
with all the required brackets, is given by: 

(Oi(ti)O 2 (t2)O3(t3)} OT0 =Tr(7r P* 1 (Oi^ 1 2 (O 2 P t t 2 3 (O3)))) • 
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D. The accompanying density matrix 

We suppose that the initial density matrix ttq of the system satisfies LqTTo = . For example, the system is 
prepared in a thermal state at temperature T = 4, where k is Boltzmann's constant and its initial density matrix 

is 7To = Zq 1 exp(—f3H(0)). If the system is closed and the Hamiltonian is constant in time, then the density matrix 
does not vary with time: 7r t = ttq- However, if the Hamiltonian H t changes with time and/or if the interactions with 
the environment are taken into account, the density matrix evolves according to equation ((2]) and the system is out 
of equilibrium. In particular, at time t, the system does not lie in the kernel of the time-dependent generator L,\. For 
example, at time t, a closed system with time-dependent Hamiltonian H t , is not in the Gibbs state Z^ 1 exp(—/3H t ). 
(The same feature is also true, of course, in classical mechanics). 

Nevertheless, given a time-dependent Lindbladian L t , it is useful, following [43], to associate to it the accompanying 
density-matrix 7r t that satisfies L\ir t = 0. Physically, 7r t represents the stationary state in a system where time is 
frozen at its instantaneous value t. For example, for a closed system, we have n t = Z^ 1 exp(— (3H t ). However, as 
we emphasized above, at time t the true density matrix of the system differs from the accompanying density-matrix: 
p t ^ 7r t . The reason is that ix t depends on time, and therefore it can not satisfy Eq. ([2]): 

— 7Tt ^ whereas L\-K t = . 



III. A QUANTUM JARZYNSKI-HATANO-SASA IDENTITY FOR LINDBLAD DYNAMICS 

A key idea that lies at the heart of the Jarzynski identity in the classical case is to consider non-equilibrium averages 
over weighted trajectories. This crucial feature was clearly recognized and stated in the very early works [13, 13 . If the 
weighting factor is chosen to be the exponential of the work performed on the system, then weighted averages along a 
non-equilibrium process between times and t can be reduced to thermal averages performed with the accompanying 
Gibbs measure at time t. An equivalent formulation due to Hummer and Szabo [45l | is to consider an auxiliary 
system governed by a fictitious dynamics, constructed in such a way that at each time t the auxiliary system lies 
in the accompanying steady-state measure at time t of the initial system. Thus, averages on the fictitious system 
can be written as steady-state averages in the initial system. Besides, using the Feynman-Kac formula, Hummer and 
Szabo showed that averages over the fictitious system are given by averages over the initial system weighted by the 
Jarzynski factor (the exponential of the work). Eliminating the auxiliary system between the two equalities leads to 
the Jarzynski identity. 

It is important to note that exactly the same idea of considering a modified dynamics appears in the proof by Kur- 
chan of the Gallavotti-Cohen relation for Langevin dynamics [54{ and also in the general derivation of the Gallavotti- 
Cohen symmetry for Markovian systems by Lebowitz and Spohn [55[ . 

Our aim in this section is to prove a Jarzynski-Hatano-Sasa identity for a quantum Markovian dynamics by using 
a similar technique: the Lindblad dynamics is deformed so that the accompanying density-matrix 7r t of the initial 
system becomes the true density-matrix at time t of the fictitious auxiliary system. Then, an operatorial version of 
the Feynman-Kac formula will allow us to relate expectations values over the auxiliary system to averages on the 
initial system, leading us to a quantum version of the Jarzynski-Hatano-Sasa identity. In the case of a closed system, 
we shall show that this identity is equivalent to relations that were previously known. Finally, we deduce from our 
general result a fluctuation-dissipation theorem valid for an arbitrary steady state. 

A. A modified dynamics for the accompanying density matrix 

The accompanying density-matrix 7r t does not obey the Lindblad equation Eq. @. However, 7r t is a tautological 
solution of the modified evolution equation 

d t n t = (i t + 7 r t - 1 (9 f ^)) t 7T t . (9) 

We introduce the non-stationary operator 

W t = -(Trt)- 1 (dtn t ) , (10) 
which reduces in the classical limit to the rate of injected power, and we define a modified superoperator as follows 

L t .i =Lt + ( 7 r t )- 1 (^t) =Lt- W t , (11) 
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where Wt acts on a given density-matrix by a multiplication on the left. Then, equation (|12p can be rewritten as 

d t % t = L{ in - (12) 

The formal solution of this equation is given by 

*t=(l$,iY*o, (13) 



the modified evolution superoperator P s 1 being defined as 

rt \ / rt 



P l sl =cx$(^J^ L u> i dv^j = exp* ^ (L u + n u 1 (d u ir u )) duj . 



(14) 



We observe that equations (fT5I [T4")) for the accompanying matrix are similar to equations ([5j for the 'true' density 
matrix of the system. 

Now consider an arbitrary observable A. Then, using equations (fl"3|). we can write 



Tr^oP^A) =Tr(7r t A) . (15) 

This identity means that averages for the fictitious evolution which are performed by using the modified evolution 
superoperator P* 1 reduce to averages performed with the accompanying density-matrix at time t. 

B. Proof of the Quantum Jarzynski-Hatano-Sasa Relation 

A quantum version of the Jarzynski-Hatano-Sasa Relation will be obtained by relating the auxiliary evolution 
superoperator P* 1 to the initial evolution superoperator P*. To achieve this aim, we need to prove an extension of 
the Feynman-Kac formula. 

We write a Dyson-Schwinger series expansion for Pq 1 , considering Wt to be a perturbation of the Lindbladian L t : 



» n 

p li = E^ 1 )" / ]> Po'w^pfrwt, . ■-p t t :_ 1 w tN p i 

„ JO<t 1 <t 2 <...<t n <t i=1 



t . 

t N ; 



(16) 



we recall that the superoperator P t ** +1 operates on all the terms to its right. This well-known formula [63[ can be 
proved by showing that both sides of the equation satisfy the same differential equation and are identical at t = 0. 
Inserting this expansion on the r.h.s. of equation (|15p . we obtain 

Tr (noP^A) = £(-l) n / J dU Tr Up^W tl P^W u _ . . . P^W^P^a) . (17) 

„ Jo<t 1 <t 2 <...<t n <t i=1 v 7 

Rewriting the trace term inside the integrals as a multi-time correlation function via Eq. (j8]), we find 

Tr^oP^A) =^(-1)" / Jdti (W tl (t 1 )W t2 (t 2 )...W tN (t N )A(t)) no . (18) 

n Jo<ti<t 2 <...<t n <t i=1 

We note that on the r.h.s. of Eq. ©, the operators Oi depend on the time tf, this time dependence is written as 
an argument Oi(U). Here we have Oi = Wti which already depends on time. A supplementary time dependence is 
introduced through the use of Eq. (|8|), which now appears as Wt^ti). By linearity, equation ([PS]) is identical to 



Tri^oP^A) = ( { / f[dU W tl (h)W t2 (t2)...W tN (t N ) \ A{t) 

\ { n Ja<t 1 <t 2 <...<t n <t i=1 J ! 

The term between curly brackets can be resummed as a time-ordered exponential 



(19) 

7TQ 



Tr (vroPoVA) = (e$ (- J W u (u) du 
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This formula is an extension of the Feynman-Kac formula for quantum Markov semi-groups. We emphasize that , due 
to non-commutativity of operators, the exponential that appears in the usual Feynman-Kac formula is replaced here 
by a time-ordered exponential. Moreover, we remark that, although there exist many generalizations of the Feynman- 
Kac formulae [l|, [25j], the present one seems to be original. Finally, using equation (fl~5j) . the following relation is 
derived: 



Tr (ir t A) = f- W u {u) du^j A(t) 



(21) 



This identity is a quantum extension of the classical Jarzynski-Hatano-Sasa identity and is one of the main results of 
the present work. In particular, if we take A = 1 then Eq. ([2~Tj) becomes 

I W u (u)du]\ =1, (22) 



\ V Jo 

where we use the fact that Tr (7r t ) = 1. Note that if we interpret the mean values as classical averages and the operators 
as commuting c-numbers, then Eq. (|2"2")l reduces to the classical Jarzynski-Hatano-Sasa relation. The relation (|2"2")l can 
be interpreted as a kind of a book-keeping formula which allows us to write identities amongst correlation functions. 
We emphasize that the operator Wt is not Hermitian: this is a signature of the fact that 'Work is not an observable' 
[70j . However, for the special case of a closed system, the work-term can be written as aproduct of two observables, 
as will be shown below, allowing us to retrieve previously known work identities @, [H, US llS [zS [Z6| ■ Besides, from 
a first order expansion of Eq. (|2"2"j) wc shall derive a generalized fluctuation-dissipation theorem valid in the vicinity 
of a quantum non-equilibrium steady state. 



C. The case of a closed quantum system 

We consider the special case of a closed, isolated system, governed by a time-dependent Hamiltonian. The Lind- 
bladian reduces to the Liouvillc operator, L t .X — i[H t ,X], and the evolution of the system is unitary. For a closed 
system, the evolution supcropcrator Pq acts on an observable X as follows 

P*X=(U t ) f XU*, (23) 
where the unitary evolution operator is defined as 

Ul = exp f du {-iH u ) . (24) 
Jo 

Here, the arrow pointing towards the left over the time-ordered exponential indicates that early times are written on 
the right and later times on the left. 

The image of X by the evolution superoperator Pq defines the Hciscnbcrg operator X H (t) where the upper-script 
H stands for Hciscnbcrg: 

P t X = X n (t). (25) 
We also note that Pq is a multiplicative superoperator i.e. for any two observables X and Y we have 

p*(xy) = p*(x)p*(y). (26) 

(This is not true in general for an open system.) Thanks to this property, the r.h.s of the general expression (J5J for 
multi-time correlations can be evaluated and we obtain: 

(0 1 (t 1 )0 2 (t 2 )...0 N (t N )) 7ro =Tr(ir O?(t 1 )O 2 K (t 2 )...O%(t N )) . (27) 

If we substitute this expression in equation ([T5]) and retrace the steps from equation (fT9|) to equation (|21[) . we find 
that for a closed system the quantum Jarzynski-Hatano-Sasa relation can be written as 



Tr (n t A) = Tr ^k q c^> (- W u (u) h du^j A H ( 



(t) , (28) 
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where from equations ([2"51 and ([2"6")l . we have 

W u {u) H = - (it- 1 )* (u) (d u ir u f (u) . (29) 

We also recall that for a closed system, the accompanying density is given by 7r f = Z^ 1 exp(—(3Ht). Equation 
can be simplified thanks to the following operator identity, which is proved in Appendix [X] 



(30) 



Substituting equation ([3"0")) in equation ([25)1 . the following (tautological) identity is obtained 

Tr (n t A) = Tr {n ((tto)" 1 ^^)) A n (t)} . (31) 

Taking A = 1, we end up with the quantum Jarzynski relation for closed systems as first found by Kurchan and 
Tasaki [H|7|: 

Tr (ttq exp(/?i7 w (0)) cxp (-/?#«(*))) - . (32) 

Hence, for closed systems, the quantum work IF characterizes a process where the energy is measured twice, at time 
and time t [70[. 

Remark: If wc suppose that H u commutes with d u H Ul then the formula (|28p can be brought into a simpler form. 
From equation ([2T))) . we have W™(u) = —d u QnZ u ) + /? (d u H u ) (u) where we explicitly used that [H u ,d u H u ] = 0. 
We then obtain 



^-Tr (ir t A) = Tr L c5$ (-0 J {d u H u ) n («) du\ A n (t)^j . (33) 
The case A = 1 gives the Hanggi-Talkner form [xO] of the quantum Jarzynski relation for closed systems: 

Tr (n cJcp- (-p J (d u H u ) n (u) d «) ) = t| ■ (34) 
Remark that this relation is true if [i7„, <9 u -ff u ] ^ 0. This fact was neglected in [7(| and corrected in an addendum of 



D. Steady-State Quantum Fluctuation Dissipation Theorem 

We now return to the general case of an open system. One main advantage of the identity (|2ip is that it implies a 
generalized fluctuation-dissipation theorem valid in the vicinity of a quantum non-equilibrium steady state. We start 
with a Lindbladian Lq, which does not depend on time, with invariant density-matrix given by ttq. We then consider 
a perturbation of Lq of the form 

L t =L - h a {t)M a . (35) 

The time-dependent perturbations h a (t) are supposed to be small and a summation over the repeated index a is 
understood. At first order, the accompanying density-matrix 7r t , that satisfies L\.ir t = 0, is given by 

ir t = 7T + h a (t)e a , (36) 

where e a satisfies 

L$.e a = Mt.Tro . (37) 

The operator W t , defined in ([TP)), reads 

W t = -h a (t)D a with D =7To 1 e„. (38) 
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We now take the functional derivative of the identity ([2"Tj) w.r.t. h a {u) for u < t. The derivative of the l.h.s. 
vanishes because irt depends only on h a (t) and not on h a (u) for u < t. We evaluate the derivative of the r.h.s. by 
using the first order expression ((38|) for W u . This yields 



6(A(t)) 



Sh a (u) 



= ±{D a (u)A(t)) vo . (39) 
h=o au 



We emphasize that the expectation value on the r.h.s. is taken with respect to the unperturbed density matrix ttq. By 
choosing At = Df,(T), Eq. (|39[) becomes structurally similar to the usual equilibrium fluctuation dissipation theorem. 
This generalizes to the quantum case a result obtained recently for classical systems [l^, |64| (see [75| for an alternative 
approach) . 

Remark: In the case of a closed system perturbed near equilibrium the steady-state quantum Fluctuation Dissipation 
Theorem (f3"9")) reduces to the celebrated Callen-Welton-Kubo relation [l2|, [52| . The details of the proof are given in 
Appendix [B] 



IV. TIME REVERSAL AND A QUANTUM TASAKI-CROOKS RELATION 

Symmetry by time reversal lies at the heart of many exact identities in non-equilibrium statistical mechanics: 
indeed, various fluctuation relations can be derived by comparing the averages performed on a given process with 
those performed on the time-reversed process. We shall first recall how time reversal can be defined for quantum 
Markov processes. Then we shall use this operation to derive a general version of the Tasaki-Crooks Relation, valid 
for open quantum systems. 



A. Time-reversal for Lindblad dynamics 



In Quantum Mechanics, time reversal on the states ip of the Hilbert space is implemented by an anti-unitary 
operator [68[ that we denote by 6. The operator 8 is anti-linear and it satisfies 9 2 = 1, 9 = 9~ 1 = 6^. For a spin-0 
particle without magnetic field, 9 can be identified with the complex conjugation operator (i.e. by time reversal, the 
Schrodinger wave function ip becomes ip*). In presence of a magnetic field, this time- inversion operation must be 
supplemented by the requirement that the reversed system evolves with potential vector A R = —A. Time reversal 
of observables (which are operators acting on the Hilbert space of states) is implemented by a superoperator K that 
acts on an operator X as follows 0, [6(| : 

KX = 9X9- 1 . (40) 

The superoperator A is anti-unitary, with A 2 = 1, K = A' -1 = and is multiplicative i.e. K(XY) = K(X)K(Y). 

We can now define time reversal for a quantum Markov process. We consider the case of a system with a constant 
Lindbladian L that lies in a steady state with density-matrix w. Conditions for defining the time-reversed quantum 
Markov process have been stated by various authors 0, H3, [H, [6(j ■ The superoperator L R that governs the reversed 
process is given by 

L R = Kn^tfirK, (41) 

where ir^ 1 and tt are understood as left-multiplication superoperators, i.e., if X is an operator we have L R (X) = 
K(tt~ 1 L^ ! (tt(KX))). The condition of micro-reversibility or detailed balance is then expressed as L R = L which is 
equivalent to Ktt~ 1 IJ'kK = L. Note also that this relation and its conjugate imply that irK = Kir, so K(ir) = tt; 
thus, the detailed balance condition also takes the form L = n^ 1 KL' ! KTr which is identical to the condition given in 

[69^ (relation 4.8). The exponentiation of this last formula leads us to the finite time formula ttPq = K (Pq) Kit 
which can be written as 

Tr (B^ttPqA) = Tr ((KA*) ttP t (KB)) (42) 

for two arbitrary observables A and B. Formula (I42[) coincides with the definition of a detailed balance given by 
Majewski [6(| inspired by Agarwal (relation 2.19 of 0) and by Fagnola et al. [38|. Besides, in [38[ a characterization 
is given of the Lindbladians that satisfy detailed balance. Finally, we must underline that there exists still another 
definition of quantum detailed balance in the sense of Frigerio et al. [M HI and Alicki which can be written in 
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the form L — 7r _1 _L^7r = 2i [H, .}. It is shown in [38| that this later definition is equivalent to the previous ones when 
the Hamiltonian and the Kraus operator are even observables, i.e. K (H) = H and K(Vi) — Vi- 

More generally, the fact that L R defines a bona fide quantum dynamics is a non-trivial fact that imposes stringent 
conditions on the initial Lindbladian L [3TL |38| . One can readily verify that the stationary density matrix associated 
with the reversed dynamics is given by tt r = Kir, because (L y(Kir) = (using the fact that LI = 0). 

Finally, we consider a non-stationary set-up with a time-dependent process, governed by a Lindbladian L t and study 
the process between the initial time t = and a final time T. We wish to consider a reversed process that also runs 
for times between and T. We emphasize that there is not a unique manner to define time-inversion, as was already 
realized in the case of classical systems [l8[ . We shall write the time-reversed dynamics by analogy with equation (|4"T]) 
and by the requirement that the accompanying distribution of the time-reversed system is the time-reversed of the 
accompanying distribution of the original system. These two conditions lead to the following Lindbladian: 

Lf, = KTX^L\-K t K with t* = T-t. (43) 

Here again, n t and 7r t _1 denote left-multiplication superoperators. Using Eq. (j43]) . and the relation L t l = we find 
(L R yKirt = 0. We thus obtain tt r = Kir t , relating, as desired, the accompanying distribution of the time-reversed 
system with that of the original system. By applying Eqs. © and © to the time-reversed system, we find that the 

corresponding evolution supcroperator of the time-reversed system is given by P*> R = ^>(^J t s du L R ^j and that the 

multi-time correlations are: 

(Oi(*i)0 2 (t 2 ) • • • N (t N )f = Tr (ifiFfr^O^^Ch ■ ■ • P^On) ■ (44) 

(The superscript R on the l.h.s. recalls that correlations are taken for the time-reversed process.) We again emphasize 
that the superoperator L R , given in equation (j43[) . must be a well-defined Lindbladian: this a non-trivial property 
This property can be ensured by imposing at each time t the quantum instantaneous detailed balance condition 
L R = L t . Here, we do not assume detailed balance and we only require the weaker condition that L R is a Lindbladian. 

B. Proof of the Quantum Tasaki-Crooks Relation 

Given a scalar a, with < a < 1, we introduce two a-deformed superoperators L t (a) and L R (a), that act on an 
observable X as follows: 

L t (a)X = (Lt + air-'dtir^X 
and L R (a)X = (Lf + aiir?)' 1 d t ir?) X . (45) 

The superoperators L t (a) interpolate between L t (the original Lindbladian) and L ti i [defined in equation (|11[)] when 
a varies from to 1. Similarly, L R (a) is an interpolation from L R to L R 1 . 
The corresponding a-deformcd evolution superoperators are given by 



Pt{a)=e$(J^duL u {a) S j and P*' R (a) = exp* (J du L R (a) 



(46) 



These modified superoperators satisfy the following key duality relation, that lies at the heart of the proof of 
the quantum fluctuation theorem (for the classical analog in which the dynamics is also modified with respect to a 
continuous parameter see [H, [HI): 

7r P T (a) - \ir T KP£' R (l - a)icj ^ . (47) 

This relation is proved by the differential equation technique. The operator Ut — kqPq (a)^ 1 is equal to 1 at t = 
and it satisfies the following evolution equation 

d t U t = TroP'ia)^ 1 (irtLtia)^ 1 - dm n^ 1 ) = U t (^L^ 1 + (<*- l)d t ir t ^) = U t (KL R (1 - a)K) f , (48) 
where we have used equation (|43[) for the last equality. Hence, we can write 



f7 r = exp*^ du(KL$_ u Ky(l-oi^J = ^K&$ dvL R {l-a)K^ = (kPq' r (1 - a)K)\ (49) 
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which proves equation (|47[) . 

Applying the duality identity (|47[) to two arbitrary observablcs A and B and using the fact that K is multiplicative 
and anti-unitary, we obtain the relation 

Tr (B\ pT(a)A) = Tr ((KA*) n R P^ R (l - a) (KB)) . (50) 

This equation is the essence of the quantum fluctuation theorem and it expresses a generalized detailed balance 
condition. [Indeed, for a system in a reversible stationary state i.e. 7r t = ir and Pq' R = Pq , it becomes identical to 
the detailed balance condition (|4"2"j) 0, [3?i H(| •] Equation ([50)1 can be brought into the following more familiar 
form by using the same method as in section UlIBI (see Appendix [Cl for more details): 

((TroBn^y (0)e5Ep- (-a /„ T du W u {u)) A(T)} = (51) 

^(tt? (KA) (Tr*)" 1 ) 1 (0)eS3 (-(1 - a)/ T duW-«(u)) (-ff-B) (T) 

where W/* denotes the reversed 'injected-power' operator, given by W/* = —(n R )~ 1 dtiT R . This identity is original 
and it implies all the other results described in the present work. In particular, if we take B = l and a = 1, then 
Eq. (|51[) is the quantum Jarzynski-Hatano-Sasa identity Eq. (|21|) . If we interpret the mean values in Eq. (|5ip as 
classical averages and the operators as commuting c-numbers, then Eq. (f5"Tj) becomes Crooks'relation. 

Remark: For a closed system, with an evolution operator Uq given in equation (|24j) . we can verify that the time- 
reversed system (|43| is also closed with Hamiltonian H$ = KH t and evolution operator Uq' R = K. (Uq) . Then, 
using identity (|4"T|) for a = 1 and the fact that K is multiplicative, we obtain 

P T (1)A = ^ 1 K [P^ R ) f K%tA = tTq 1 KUq ,r K (tt t A) (u^ R ) ' = (52) 
n^K (K ((^) f ) K (n T A) K (U*)) = ^ ([/ r ) t 7T T AC/ r . 
Substituting the last expression in Eq. ([50]) leads to 

Tr (Bt^^-i^t^ T A;7 T ) = Tr (V (^) ^ (u^ R ) ^ K (B) . (53) 

Recalling that tt^ 1 and 7Tt arc given by the Boltzmann law, the above equation becomes in the Heisenberg represen- 
tation denoted by the superscript H, 

Tr{B^ o exp(pH^(0))exp(^m¥(T))A H (T)) = ^-Tr (k (rf) k r K (B) . (54) 

We emphasize that for B = 1, Eq (|53[) is a tautology (because K is anti- unitary ) , however it implies the non-trivial 
result ([54]) : this feature is characteristic of most of the derivations of the work identities. If we take A = B = 1, we 
retrieve the quantum Jarzynski relation for closed systems as first found by Kurchan and Tasaki (53l . [73| . 

V. CONCLUSION 

In this work, we have derived fluctuations relations for an open quantum system described by a Lindblad dynamics 
that takes into account the interactions with the environment as well as measurement processes. We prove the 
fluctuations relations thanks to a suitable deformation of the system's dynamics: this crucial technical idea provides 
a truly unified picture of the fluctuations relations, whether classical or quantum, and does not require to define 
the concept of work at the quantum level. Quantum Fluctuation Relations for open systems are, at present, not 
as developed as their classical counterparts. One major difficulty in the quantum realm is the lack of a trajectory 
picture when coherence and measurements are taken into account. To overcome this difficulty, an unravelling of the 
Lindblad equation [HIHH must be used. Previous attempts of this idea were performed in (26|, [29|, [3(| and we plan to 
extend the results of the present work by using such unravellings @ . Another possible extension of our work would 
be to study the effect of choosing a time inversion different from that of Eq. (|43| : this may lead to various families of 
fluctuation relations, as happened in the classical case [l8j . One could also study particles with non-zero spins such 
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as Dirac spinors. Exact solutions of specific models (such as quantum Brownian motion) may also provide us with 
experimentally testable predictions; more precisely, the formal relation between the classical exclusion process and 
the Lindblad evolution of free fermions in one dimension [32l . l74l |. could allow us to use, for open quantum systems, 
the exact results obtained for the large deviation functions of various stochastic processes [13, UK]. Finally, the 
investigation of time-reversal properties of quantum Non-Markovian systems, in which the characteristic time scale 
of the environment can not be neglected w.r.t. that of the system, should also yield interesting fluctuation relations. 

R. C. thanks K. Gaw§dzki for pointing out the fact that the Lindbladian character of Eq. (|43)l is non-trivial and 
the relation with detailed balance. R. C. acknowledges the support of the Koshland center for basic research. K.M. 
thanks M. Bauer and H. Orland for useful comments and S. Mallick for useful remarks on the manuscript. Results 
similar to those presented here were also reached independently by K. Gawgdzki and S. Attal some time ago 



Appendix A: Proof of Equation (|30|) 

The identity (|30l) is proved by using the differential equation technique. First, we note that both sides of equa- 
tion (|3T))) coincide at t — 0. Then, we find that the time derivative of the l.h.s. is given by 



d_ 

dt 



l&$ (- J W u {u) H dM = (- J W u ( 



(u) H du\(-W t (t) H ) 



This follows from the very definition of a time-ordered exponential. 
The time derivative of the r.h.s. is given by 



= hr^J^t))" 1 (if)"} = (n Q )-\?(t)(-W t (t)K) (Al) 



where we have used the fact that 7r t commutes with H t . The last equality follows from the definition ([10]) of Wt ■ We 
have thus shown that the l.h.s. and the r.h.s coincide at t = and that they satisfy the same first order differential 
with respect to time: they are therefore identical for all times. 



Appendix B: Proof of the Callen-Welton-Kubo Formula using the steady-state quantum fluctuation 

dissipation theorem ()39l) 



In this Appendix, we show that the steady-state quantum Fluctuation Dissipation Theorem (1391) is equivalent, in 
the case of a closed system perturbed near equilibrium, to the Callen-Welton-Kubo relation [l~2L l52|. 

We start with a time-independent Hamiltonian Hq with invariant density-matrix given by ttq = Z^ 1 cxp(— (3Hq). 
We then consider a perturbation of H of the form 

H t =H Q -h a (t)O a . (Bl) 

In a closed system, the Lindblad equation reduces to the Liouville equation 

d tPt =L\ Pt = -i[H t ,p t }. (B2) 

The accompanying density is explicitly given by 7r f = Z^ 1 cxp(— (3H t ). Comparing with equation (|35[) . we find M a = 
i[O a ,.]. 

We now derive an explicit expression for the operators D a defined in equation (|38p . Starting with the exact formula 
exp(-PH t ) =ex$\h% daexp(-aiJ)O a exp(ai/)^ exp(-/?if) , (B3) 
(which can be proved by differentiating both sides w.r.t. /?), we find at first order (63j | 

rP 

exp(-pH t ) = exp(-/3H) + h a t da exp(-a#)O a exp(aiJ) exp(-/3£T) + o{h) . (B4) 

Jq 
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This implies that 



Z t = Tr(exp(-pH)) + h%/3Tr(O a exp(-/3H)) + o(h) 

= Z (1 + hW (Oa)J + 0{h) . 



(B5) 
(B6) 



The first order perturbation of Wt is then given by 



W t = h1\P{O a ) -exp(+/3i?) / daexp{-aH)O a exp(aH)exp(-f3H) ) + o{h) 



Comparing with equation (|38[) . we obtain the analytical expression for D a : 



D a = -p (O a ) + f da cxp(aH)O a exp(-aH) . 
Jo 



(B7) 



We now transform, using ©, the r.h.s. of the quantum Fluctuation Dissipation Theorem (|31?)) as follows: 

£ (D a (u)A(t))„ o = A-Tr(n D a P^A) = -Tr(n D a LP^A) 
= -{Dl-K ,LP t u A) = -{tf(Dl-Ko),P t u A) 
= -Tr((Lt (D^o)) 1 PlA) = -Tr^ 1 (tf (n D a )) I* A) . 

Note that in the first equality, we use the fact that 7To is the invariant density of the unperturbed dynamics. The 
quantum Fluctuation Dissipation Theorem (|39p can thus be rewritten as 



S(A(t)) 



(B8) 



h=0 



(B9) 



Sh a (u) 

where we have defined 

E a = -TTq 1 ^ {lT D a ) 

From equation (|B7[) . we deduce the analytical expression of E a : 
E a = -ttq 1 ^ {n Q D a ) 

f P 

= icxp((3H) H,exp(-PH) / dacxp(aH)O a cxp(-aH) 

Jo 

rP rP 
H, / da exp(aH)O a exp(— aH) =i daexp(aH)[H,O a ]exp(—aH) 
Jo Jo 

r p d 

= % I da— (cxp(aH)Oa cxp(-aH)) = icxp(/3H)O a exp(— (3H) — iO a . 
Jo da 

We remark that the terms on the r.h.s. can be interpreted as the analytic continuation in imaginary time (as allowed 
by the KMS condition ;59|) of the Hcisenberg representation with respect to the unperturbed Hamiltonian H. Thus, 
we have 



E a = iOf{-iP)-iOf{Q) 



(BIO) 



Finally, the (|B8| becomes for (it < T) 
5{A{t)) 



Sh a (u) 



h=0 



= iTr (Tr O^(-iP)A H (t - u)) - iTr (ttqO* (Q)A h (t - u)) 

= i (O a (-i{3)A(t - u)) Wo - i (O a (0)A(t - u)) no 
= i (O a (0)A(t -u + 03))^ - i (O a (0)A(t - u)) no 



(Bll) 



14 



The last equality follows from the fact that the correlation (X a Yt) depends just on t — s. This equation is the real 
space version of the Callen-Welton-Kubo equation. The more conventional form is obtained by performing a Fourier 



Transform with respect to time. The susceptibility is defined as (using causality i.e. -gj^rr 

8{A(t)) 



h=0 



if u > t) 



Xo a A{w) 



(It 



8h a (0) 



cxp(iwt) 



h=0 



i(exp((3w) - 1) / dt (O a (0) A(t)) exp(iwt) 



(B12) 



where we have used equation (jBlip in the last equality. The symmetrized correlation can be written as 



C S 0aA {w) 



dt exp(iwt) Tr(n OZ (0)A H (t)) 



1 



— OO 



dtcxp(iwt)Tr(-K OZ (0)A H (f)) + 



1 

2 7, 

1 



— oo 



dtcxp(iwt)Tr(ir A H (t)O^ (0)) 
dtexp(-iwt)Tr{n O^ {0)A H (-tj) + 



+ 00 



dtexp(iwt)Tr(Tr A H (t)O^ (0)) + - / dt cxp(-iwt)Tr(ir A H (-i)Of (0)) 



+ 



dtexp(iwt)Tr(-K O" (0)A U {t)) + - / dt exp{-iwt)Tr{-K A u (0)O„ (* + ifi)) + 



dtexp(iwt)Tr(-K OZ (Q)A H (t + i0)) + 



+ 'X 



2 7o 



dtexp(-iwt)Tr(-n A H (O)O^ (t)) 



= -(l + exp(pw)) di(O„(0)A(t)) exp(iu;i) + -(1 + cxp(-/3w)) / dt (A(0)O a (t)) exp(-iurf)(B13) 
2 Jo 2 

where we have used a KMS type identity to obtain the third equality: Tr{it Q X H (s)Y H (t)) = Tr{it Q Y H ~(t)X H '(s + i/3)). 
By using (|B12j) and (|B13|) . we obtain 



Cm " 1 = ^(l + exp(/3w t-— — - + - 1 + cxp (-f3w — — — — 

2 z(exp(/3u;) — 1) 2 ?(exp(— pw) — 1) 



(B14) 



coth(^) 
2i 



(XO a A(w) - XAO a (-w)) 



Note that in the expression for XAO a , the operator A is considered to be the perturbation and O a the observable. We 
thus implicitly suppose that A = A'. In particular, if we take A = Ob 



c s o a0b H 



coth(^) 

(xo a O b (w) - Xo„o a (-W)) = coth(— ) Im(xo a O b (w)) . 



2 1 



~2~ 



(B15) 



recalling that XAO a {~ w) = XAO a {w). For a closed system, an alternative proof showing that relation (|3ip implies 
Eq. (|B15|) is given in Q. 



Appendix C: Derivation of Equation (|51[) 



First, we establish the following relation, valid for two operators X and Y 

Tr(n YpT( a )X) = (y(0) cx^ f-ajf dKW u (u)j X(T)\ . (CI) 

The method to derive this formula is identical to the one used in equations (|16[) to ((20]) : we perform Dyson-Schwinger 
expansion Pq(o.) w.r.t. the deformation parameter a, rewrite the trace as a correlation function via Eq. (JSJ and the 
result is rcsummcd as a time-ordered exponential. We recall that Wt was defined in equation (flU)) . The same equation 
is true for the reversed system, with a replaced by 1 — a and where the reversed 'injected-power' operator is given by 
W t R =-{ir*)- 1 d t *?: 



Tr (n YP^' T (l - a)x) = (y(0) cx^ ^-(1 - a) du W* (u)j 



X(T)) . (C2) 
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Replacing both sides of equation ([50)1 by the relations (jCll) and (JC2J) and inserting the duality identity (jT?)) into 
equation f|Cl[) . leads to the Fluctuation Relation ([5"Tj) for an open quantum Markovian system. 
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